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ABSTRACT 

We present single-epoch multi-wavelength optical-NIR observations of the "anoma- 
lous" lensed quasar SDSS J0924-I-0219, made using the Magellan 6.5-metre Baade 
telescope at Las Campanas Observatory, Chile. The data clearly resolve the anoma- 
lous bright image pair in the lensed system, and exhibit a strong decrease in the 
anomalous flux ratio with decreasing wavelength. This is interpreted as a result of 
microlensing of a source of decreasing size in the core of the lensed quasar. We model 
the radius of the continuum emission region, cr, as a power-law in wavelength, a cx 
We place an upper limit on the Gaussian radius of the u'-band emission region of 

cm, and constrain the size-wavelength power-law index 
to C < 1.34 at 95% confidence. These observations rule out an alpha-disc prescription 
for the accretion disc in SDSS J0924-I-0219 with 94% confidence. 

Key words: accretion discs gravitational lensing quasars: individual: 
SDSS J0924-F0219 



1 INTRODUCTION 

Half a century on from the discovery of quasars our 
paradigm for understanding their central engine remains 
poorly physically motivated and only very weakly observa- 
tionally constrained (|Blaes'2007!). While the standard accre- 
tion model has met significant success in describing Cata- 
clysmic Variables and Young Stellar Objects, it does a poor 
job when dealing with accretion onto black holes, and is ther- 
mally and viscously unstable in situations where radiation 



pressure dominates gas pressure (see the review by Blaes 
|2004 ). AGN accretion discs are impossible to observe di- 
rectly, having typical angular sizes on the order of micro-to- 
nano arcseconds. 

Microlensing is an established technique for obtain- 
ing high spatial resolution measurements of lensed objects 



(e.g. Wambsganss and Paczynski 1991| Rauch and Bland- 



fordl|1991| [Witt et al.||1995| |Wyithe et al.||2000[ [Schechter 



and Wambsganss 2002' ), but until recently the technique 



was not refined enough to place meaningful statistical con- 



straints on lensed quasar emission regions. Bate et al. ( 2007 1 



presented a technique for calculating source size probabili- 
ties using Monte Carlo simulations of magnification maps 
at the minimum and saddle-point image of an "anomalous" 



* This paper includes data gathered with the 6.5 meter Magellan 
Telescopes located at Las Campanas Observatory, Chile, 
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(NFB), rwebster@unimelb.edu.au (RLW) 



quadruply-imaged lensed quasar, which allows us to place 
far tighter constraints on the central source sizes of quasars. 

In an AGN accretion disc, gravitationally bound mate- 
rial spirals inwards onto a central black hole. Gravitational 
compression heats the material, which then radiates. In or- 
der to fall inwards, matter must lose not only energy but also 
angular momentum, while the total angular momentum of 
the system is conserved. The central problem in accretion 
discs is one of angular momentum transport outwards in 
the disc, in order that material can continue to plunge in- 



wards. See the excellent reviews by Blaes (20041 and Balbus 



(2003) for detailed information. While molecular viscosity 
helps heat up the disc and thus radiate energy, it is insuf- 
ficient to transport the angular momentum outwards. The 
idea of an anomalous stress introduced by nonlinear physics 
("turbulence") dates back to|Shakura and Sunyaev (19731. 



A longstanding problem has been how to introduce such 
a turbulence to what is essentially a stable laminar flow - 
the Rayleigh Stability criterion excludes any hydrodynamic 
mechanism (e.g. |Frank et al.|2002|rBalbus|2003]|Blaes|2004 1. 
A range of theoretical models based on different physical 
mechanisms have predicted different radial emission profiles. 

[Shakura and Sunyaev| ( |1973[ ) first invoked a non-linear 
turbulence to produce an additional stress in the disc, pa- 
rameterised by a coefficient < a < 1. Based solely on 
energetics and conservation of mass (ignoring the issue of 
angular momentum transport entirely), one can derive a 
simple radial temperature profile, Tc oc r~^, with p = 3/4 
(see e.g.|Pringle and Rees|1972[ [Shakura and Sunyaev|1973| 
Gaskell 2008 1 . If the emergent spectrum is the super-position 
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of blackbodies generated locally, T oc \rn\x (e-g- Frank et al 



yaev ( 1973 1 disc scenario in this source. Note that simi- 



20021. This gives a spectral profile 

r oc A"'^'''' oc 



(1) 



where ^ = 4/3. This is what you get from a simple Shakura 



and Sunyaev (19731 treatment of the disc. Gaskell (20081 



suggests the opposite extreme, in which all energy simply 
comes from the outward flow from hotter to cooler regions. 
With a single central source, and thus a flux that falls off as 
r~^, the effective temperature drops off as T oc £~^^'^ 
P 



giving 



0.5, oc u ^ and C — 2. The "slim disc" (Abramowicz 
et al.|1988 1 is an example of a model that gives rise to such a 



temperature proflle. Gaskell ( 2008 1 also calculates the value 



of p required in order to flt a generic observed quasar spec- 
tral index in the optical-UV: F^, oc f""'^. This requires a 
value of p = 0.57, and thus — 1.75. It has also been sug- 
gested that magneto-rotational instability (MRI) might be 



a generic source of turbulence in discs ( Balbus and Hawley 



1991 1. Agol and Krolik ( 2000 1 present an example of such a 



model that assumes a strong magnetic connection between 
black hole and disc that essentially "spins up" the disc from 
the black hole's own rotation. This gives rise to a temper- 
ature proflle with p = 7/8, or a wavelength proflle index, 
C = 8/7. 

Anomalous lensed quasars have an image pair in 
which one of the images is unusually ( "anomalously" ) dim, 
with the evidence strongly suggesting that the anomaly is 
caused by microlensing. Observations have shown that the 
brightness anomaly increases with decreasing wavelength 



(e.g. Wambsganss and Paczynski 1991 Inada et al. 2003 



Bate et al.||2008| |. Some studies of such sources have at- 



tributed the anomaly to dust, e.g. Lawrence et al. (19951; 



McLeod et al. (19981; 0stman et al. (20081. Most of these 



studies were carried out before the relevant lensing galaxy 
redshifts were known, and more recent observations gener- 
ally cast doubt on the dust hypothesis. In most (eight out 
of ten) known cases ( [Pooley et al.|[2"007[ ) the anomalously 
dim image lies at the saddle-point of the Fermat light-travel 
time surface as predicted by lensing theory. It is difficult 
to achieve this with dust (P^g/io = 0.04 from the bino- 
mial distribution). Furthermore, X-ray observations of the 
quadruply lensed quasar IRXS J1131-1231 show a very low 
absorbing column density, ruling out differential absorption 
in at least this case ( Blackburne et al.|2006 1 . Millilensing by 
galactic substructure has also been explored, but would not 



produce the observed polychromatic effect (IMao and Schnel 



derl[T998l [Metcaff and Madau|[200T| |Chiba]|2002| [Dalai and' 



Kochanek|2002[ ). Of course gravitational lensing is a funda- 
mentally achromatic process at any scale, and the microlens- 
ing interpretation introduces a chromatic effect by virtue of 
the lensed source size varying with wavelength. A strongly 
anomalous radio flux ratio would provide strong evidence 
for millilensing, but none has been observed to date. 

Assuming that the anomaly in our flrst observed ob- 
ject (MG 0414-1-0534) is entirely due to microlensing, |Bate| 



et al. ( 2008 \ placed some of the strongest constraints yet 



on accretion disc size, but placed only weak constraints on 
the emission mechanism, excluding the simple turbulence- 



driven Shakura and Sunyaev ( 1973 1 model at only the 68% 



level. Here we present our second case, SDSS J0924-I-0219. 
We are able to place much stronger constraints on the ob- 
served index and effectively rule out the [Shakura and Sun-| 



lar work can be done using long-term monitoring of non- 
anomalous lensed quasars ( Eigenbrod et al^|2008 L The ad- 
vantage of using anomalous lensed quasars is that the close 
proximity of the anomalous image pair guarantees a short 
relative time delay between the two images. Thus observing 
the sources close together in time means we are not seeing 
the source at two different epochs. Our technique is appli- 
cable to only a small subset of known lenses, but for this 
subset, using single epoch data, tight constraints can be ob- 
tained, that in other sources would take decades. 

We assume throughout this paper a cosmology with 
Ho = 70 km s"^ Mpc"\ r2,„ = 0.3 and = 0.7. 



2 OBSERVATIONS 

SDSS J0924-I-0219 was the second lensed quasar to be dis- 
covered by the Sloan Digital Sky Survey (SDSS). It is the 
most anomalous lensed quasar known. It was previously im- 

(|2003|) in 0'.'55- 



Inada et al. 



aged on Magellan (MagIC) by 
0'.'75 seeing as a follow-up to initial discovery by SDSS. The 
anomalous pair of images, labelled A (brightest image) and 
D, are weakly resolved in these data. We obtained the orig- 
inal Magic imaging data taken by [Inada et al.[ ( |2003| ) in 
order to work directly on the raw images. There are no ex- 
isting radio, mm nor sub-mm detections of the source, and 
no mentions of dust in this source in the literature. IKeetonI 



et al.[ ( |2006| ) presented HST/ACS and NICMOS imaging 
and spectroscopy of the source, showing that the anomaly is 
present in both the continuum and the broad emission line 
flux ratios. Chandra soft X-ray measurements also show a 

The redshift of the 



20071 



strong anomaly (Pooley et al. 

quasar is Zs = 1.524 ( [Inada et ia,l.[[2063| and the red shift of 
the lens galaxy is zi = 0.394 ( Eigenbrod et al.|2006 |. 



In this paper we present a third epoch of optical-NIR 
observations, taken with the Magellan 6. 5- metre Baade tele- 
scope on 2008 March 21 in 0'.'3 - 0'.'45 seeing. These ob- 
servations complement the earlier two epochs of optical-IR 
data, spanning a larger wavelength range (0.4 - 1.6/im) near- 
simultaneously, and in greater spectral resolution (seven 
bands) than any earlier dataset. Observations were made 
using IMACS (in f/2 mode) in SDSS g',r',i' and z' bands, 
and PANIC mY,J and H. Integration times were approx- 
imately flve minutes per fllter. Images are provided in Fig- 
ure [l] Photometric calibrations were performed using three 
standard star observations over the course of the night in 
the optical, and two in the NIR. We note that any biases in 
the calibration will not affect the D/A ratios, and thus our 
modelling of the data, which depends only on the relative 
fluxes of the anomalous image pair A and D. Wc fit Gaus- 
sians to each quasar image, beginning with an assumption 
of a round, Gaussian PSF, and then allow the ellipticity to 
vary as a free parameter. We initially used the HST data 
to constrain the position of each component in each of our 
observations. In the final fit we freed up position as well, 
resulting in a 5-D fit (x, y, amplitude, FWHM, and elliptic- 
ity) for each component using a downhill gradient technique. 
The apparent magnitudes of the four lensed images and the 
lensing galaxy are provided in Table [T] Flux ratios in each 
filter between the two images of interest. Fobs ~ D/A, are 
provided in Table [2] This table contains flux ratios from 
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the observations presented in this paper, as weU as Chandra 
data ( Pooley et al.|[2007 l, HST observations obtained from 
the CASTLES Survey web pagaj and the earlier Magellan 
observations of Inada et al. ( 20031). 



As discussed by Keeton et al. (20061, microlensing 



3 MODELLING 

We parameterise the radius of the emission region cr as a 
power-law in wavelength in order to test against various 
accretion-disc models: 



A 



(2) 



We use the modelling technique desc ribed by |Bate et alT] 
( 2007| ), and applied to MG0414-I-0534 ( |Bate et al.|2008| ). To 
re-iterate briefly, we begin with a macro-lensing model ob- 



tained from Keeton et al. (20061. This provides the conver 



gence k and shear A at each image position. We consider two 
distinct components to convergence - a continuous compo- 
nent Kc and clumpy, compact stellar component At*. Caustic 
networks are then constructed using an inverse ray-shooting 
method ( Wambsganss et al. 1990 1 . The microlens population 



is drawn from a Salpeter mass function, dN/dM oc M~^'^^, 
with Mmax/Mmin = 50. Physical sizes are scaled by the av- 
erage microlens Einstein radius projected onto the source 
plane, rjo = 5.78 x lO^'^h'^^'^ {{M)/Mq)^^^ cm. We also per- 
formed simulations using a flxed mass microlens population 
(at the average mass of the Salpeter function above). 

Gaussian source surface brightness profiles b of various 
characteristic radii a are convolved with the caustic net- 
works in the image plane: 



b = exp 



2(j2 



(3) 



We then construct a library of simulated observations of fiux 
ratio as a function of source radius. These are converted 
to flux ratio as a function of wavelength using Equation [2] 
Comparing this library of simulations with the observational 
data allows us to build up three dimensional probability dis- 
tributions for smooth matter fraction, radius of the inner- 
most emission region ao, and power-law index ^. Smooth 
matter was allowed to vary between 0% and 99% in 10% 
increments. We obtain the results presented in this paper 
by marginalising over smooth matter fraction: 



d^P 
dC^dao 



d-^P 
dsd(dao 



ds 



(4) 



Source sizes were 0.05 770 to 2.00 t^o, in steps of 0.05 
rio- The upper end is large enough that microlensing should 
no longer be significant. The lower end is just above the 
resolution limit of our simulations, and is equivalent to the 
Schwarschild radius = 2GM/c^ of a 10^ Mq black hole. 
The mass of the black hole in SDSS J0924-I-0219 has been 



estimated at between 1.1 x 10® (iPeng et al. 20061) and 



2.8 X 1O*M0 ( [Morgan et al.|2006| ), so the lower limit of our 
calculation is close to the expected innermost stable circular 
orbit of a nonrotating black hole, ZRs — QGM/c? . 



http: / / cfa-www.harvard.edu / castles / 



alone is sufficient to explain the observed anomaly. How- 
ever, Eigenbrod et al. (2006) have argued for the presence 



of CDM substructure in the vicinity of image D. Thus 
millilensing may also be present. However, this would not 
produced the observed chromatic effect, rather it would off- 
set the total observed anomaly. 



4 RESULTS 

Our newly acquired images of SDSS J0924-I-0219 are pre- 
sented in Figure [1] In Table [l] we present our PSF photom- 
etry. Table [2] shows that the fiux ratio between images D 
and A grows larger as wavelength increases. We also show 
all published values of the anomalous fiux ratio in the liter- 
ature, from X-ray to NIR. The anomalous flux ratios are 
plotted against wavelength in Figure [2] We assume that 
this anomaly is solely due to microlensing, and apply mi- 
crolens Monte Carlo simulations as described above to re- 
turn a joint probability distribution in source size, radial 
wavelength power index and smooth matter fraction in the 
lens. Marginalising over the latter we obtain the probability 
distribution presented in Figure |3] Cumulative probability 
distributions in source size a and radial wavelength index C, 
are also shown. We model the size of the continuum emission 
region, a, as a power-law in wavelength, ct cx A"- . We place an 
upper limit on the Gaussian radius of the emitting region of 
3.04X lO'''/i7o'''^({M)/M0)i/2 cm {\1.7?.h-^^^{{M) /MqY'^ 
light days) in it'-band, and constrain the power-law index to 
C < 1.34 at 95% confidence. 

If we add in all the observational data (assuming that 
each epoch of data applies an independent constraint to the 
accretion disc structure) we obtain the slightly tighter prob- 
ability distribution shown in Figure |4] Together they con- 
strain the radius of the emitting region to be < 1.02 x 
10"/i-'/''((M)/Mq)1/2 cm (3.95fe-'''^({M)/M0)i/2 light 
days) in u'-band, and the power-law index to 0.06 < C, < 
1.24 at 95% confidence. 

We have presented results obtained using a Salpeter 
microlens mass function as we feel they are most gen- 
erally applicable. Switching to a fixed mass microlens 
population (mean mass of the Salpeter function) results 
in near-identical constraints on (^ < 1.32), but with 
slightly tighter constraints on the actual source size: O.SIt^o, 
or 1.77 X lG^^'h~^''^{{M)/MQY^'^ cm in ii'-band , or 
&mh~Q^'^{{M)/MQY'^ light days at 95%, using the Magel- 
lan data only. 



5 DISCUSSION 



Keeton et al. ( 2006 ) established that microlensing is present 



in this system, based on the difference in the continuum and 
broad emission line fiux ratios and its photometric variabil- 



ity on the timescale of a decade (Kochanek 2005 1. Eigenbrod 



et al. (2006) see rapid (15 days) and asymmetric changes 



in the emission lines of the quasar images, and no strong 
changes in the continuum. We find the time delay between 
the two close images A and D to be less than a day, using 
the adopted macrolens mo del of [Keeton et al.| ( [2006| l and 
the GRAVLENS software ([Keeton|[2001[). See also table 3 
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Table 1. New SDSS J0924+0219 photometry from Magellan IMACS/PANIC: Apparent magnitudes of the lensed images (A, D, B, and 
C) and the lensing galaxy (G) in SDSS J0924+0219. Observations were taken in seven filters with the IMACS and PANIC instruments 
on the Magellan 6.5-m Baade telescope, 2008 March 21. Typical photometric errors are 3.5% in the NIR and 6% in the optical. 



Filter 


A 


D 


B 


C 


G 


H 


17.39 ± 0.04 


18.99 ±0.04 


18.34 ±0.04 


18.80 ±0.04 


17.39 ±0.04 


J 


18.24 ±0.04 


20.29 ±0.04 


19.44 ±0.04 


20.02 ±0.04 


18.16 ±0.04 


Y 


18.30 ±0.04 


20.46 ± 0.04 


19.51 ± 0.04 


20.09 ±0.04 


19.07 ± 0.04 


z' 


19.84 ±0.06 


21.65 ±0.06 


20.98 ±0.06 


21.65 ±0.06 


21.15 ±0.06 


i' 


19.54 ±0.06 


21.54 ±0.06 


20.77 ± 0.06 


21.40 ±0.06 


21.06 ±0.06 


r' 


19.85 ±0.06 


22.32 ±0.06 


21.31 ± 0.06 


22.59 ±0.06 


20.41 ± 0.06 


g' 


20.14 ±0.06 


22.90 ±0.06 


21.57 ±0.06 


22.72 ±0.06 


21.63 ±0.06 



Table 2. Central wavelengths and observed (anomalous) flux ra- 
tios -Foi,3 between images A and D in each filter. The 2008 March 
21 observations were taken with the IMACS and PANIC instru- 
ments on the Magellan 6.5-m Baade telescope. '"'2001 December 
15 data were taken using MagIC on Baade ( |Inada et aT1|2003| . 

2003 data were t aken using HST/NIC MOS and WFPC2 for the 
CASTLES survey ("Keeton et al. '20061. ('='2005 data were taken 
using Chandra ( (Pooley et al.||2007| ). 



Filter 


Ac (A) 


Fobs = D/A 


Date 


H 


16500 ± 1450 


0.23 ±0.05 


2008 March 21 


J 


12500 ± 800 


0.15 ±0.05 


2008 March 21 


Y 


10200 ± 500 


0.14 ± 0.05 


2008 March 21 


z' 


9134 ± 800 


0.19 ±0.10 


2008 March 21 


i' 


7625 ± 650 


0.16 ±0.10 


2008 March 21 


r' 


6231 ± 650 


0.10 ±0.10 


2008 March 21 


9' 


4750 ± 750 


0.08 ±0.08 


2008 March 21 




7625 ± 650 


0.08 ±0.05 


2001 December 15 


j.'(a) 


6231 ± 650 


0.07 ± 0.05 


2001 December 15 




4750 ± 750 


0.06 ±0.05 


2001 December 15 


n't") 


3540 ± 310 


<0.09 


2001 December 15 


Fieow^'' 


15950 ± 2000 


0.08 ±0.01 


2003 November 18 


F814W(*') 


8269 ± 850 


0.05 ±0.005 


2003 November 18-19 


F555W('') 


5202 ± 600 


0.05 ±0.01 


2003 November 23 


0.5-8keV('=) 


12 ± 11 




2005 Feb 24 



high resolution radio imaging we cannot rule out millilens- 
ing entirely. A strong anomalous ratio D/A in the radio 
would strongly confirm the presence of millilensing. If con- 
firmed this would shed interesting hght on the nature of dark 
matter in the lens, while not strongly affecting our conclu- 
sions, since the observed chromaticity cannot be explained 
by millilensing. 

What happens if we assume that the anomaly is par- 
=±tially or wholly dust? There are three strong arguments 
against differential extinction by dust in this source. Firstly, 
we would expect a far smaller effect in the NIR, where little 
extinction would be expected. To get such strong NIR ex- 
tinction, we would expect extraordinarily large extinction 
in the optical. This is illustrated in Figure [2] No physi- 
cal dust model exists that can produce the observed effect. 
Dust cannot explain the shape of the observed drop off with 
wavelength. Secondly, in eight out of ten known anomalous 



lenses, it is the saddle image which is fainter (see Pooley 
|et al.|2007 l. The probability of dust extinction would be the 
same for each image. And finally, in one other case, extinc- 
tion by dust has been has been shown to be negligible from 
the x-ray observations (Blackburne et al. 20061). If we as- 



sume that millilensing accounts for all of the observed effect 
in the NIR, and that dust caused the remaining differential 
.extinction, we would still see a far steeper drop off in the 
flux of image D to the blue. 



of Eigenbrod et al. ( 2006 1 who present the time delays for a 
variety of macrolensing models, all of which are consistent 
with IaD < 1 day. Our results therefore simply rely on the 
choice of the model parameters for the particular images in 
this source. 

Unfortunately in their follow-up Keck spectra |Inada| 
~ 1^03) observed A and one of the other (non- 



et al. 



anomalous) images, B. However, [Keeton et al.| ( |2006[ ) 
present spatially resolved HST spectroscopy of the entire 
system. The authors note that there are differences between 
the emission-line and continuum flux D/A ratios. The D/A 
ratio is 0.1 in the broad Lya line but 0.05 in the associ- 
ated continuum. They argue that microlensing can account 
for both the continuum and emission-line flux ratios, if the 
broad emission line region is comparable in size to the Ein- 
stein radii of the microlenses. Comparing our flux ratios to 
those of Inada et al. and Keeton et al. conflrms that we are 
indeed seeing the effects of microlensing, since the values 
have changed significantly on a 5 year time scale. However, 
all the existing NIR ratios also show a significant anomaly. 
There are no radio detections for this source and without 



6 CONCLUSIONS 

We have modelled the size of the continuum emission re- 
gion, a, as a power-law in wavelength, cr oc A"- for the lensed 
quasar SDSS J0924+0219 . Using our new Magellan data 
alone we place an upper limit on the Gaussian radius of the 
it'-band emission region of 3.04x 10^'^/i7q^''^((M)/Mq)^/^ cm 
(a little under 12h:^Q^'^ {{M) /Mq)^^^ light days), and con- 
strain the power-law index to ^ < 1.34 at 95% con- 
fidence. This is tighter than the constraint placed on 



MG0414+0534 ( [Bate et al.||2008| ), and excludes the ac- 
cretion model of |Gaskell| ( |2008|) at greater than th e 99% 
confidence level. The simple Shakura and Sunyaev ( 1973 1 



disc is excluded at lower confidence (94%). Combining our 
data with earlier epoch HST (ACS and NICMOS) and 
Magellan (MagIC) data allow us to exclude all alpha- 
disc models ( [Shakura and Sunyaev 1973 Gaskell 2008 \ 
at 97% confidence, and constrains the source size to < 
4h^o'''^{{M)/MQy/^ light days. 
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We explore extinction of image D by dust, and find a 
number of arguments against this possibility. Tiiis hypothe- 
sis could be comprehensively tested with high resolution NIR 
spectroscopic observations aimed at resolving any metal ab- 
sorption lines. Millilensing remains a possibility that would 
explain the extremely anomalous NIR flux ratio, and which 
could be ruled out with VLBI radio observations. However 
millilensing cannot explain the observed chromaticity of the 
anomalous flux ratio, and even millilensing and dust to- 
gether do not explain the observed spectral shape of the 
source. 

Results for two sources from the anomalous lens tech- 
nique are consistent. In Figure [5] we compare published 
lensing constraints on quasar accretion disc mechanisms 
from the literature. It can be seen that there is a consis- 
tency amongst techniques that is apparently converging on 



the Agol and Krolik ( 2000 1 model 



Data from additional new sources is required to fur- 
ther strengthen the results. Only a small number of such 
serendipitous targets exist, but if the current trend continues 
for other sources, we will be able to place strong constraints 
on the radial profile of the accretion disc in quasars. 
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Figure 2. Anomalous flux ratio with wavelength for this study 
(filled squares) and for earlier data: Circles |Inada et al.| ||2003|; 
Triangles [Keeton et al^ | |2006[ l; Down-pointing triangle | |Pooley| 
|et al.|2007[ l. The dashed line shows the best quadratic fit to our 
new data only. For reference, we model the ratios if we assume 
that D is simply extinguished by dust at the redshift of the lens, 
zi = 0.394 | |Inada et~ al. 2003), with no microlensing. Two dust 
models are shown jMathis 1990) with Ry = 3.1 (dotted line - 
typical diffuse ISM) and Ry = 5 (dot-dashed line — dense dust 
clouds). Partial extinction would result in the lines moving to 
the right, with the same slopes. Millilensing can help explain the 
strongly anomalous NIR fiux ratio, but cannot explain the slope. 



Eigenbrod, A., Courbin, F., Meylan, G., Agol, E., Anguita, 
T., Schmidt, R. W., and Wambsganss, J.: 2008, A&A 490, 
933 

Frank, J., King, A., and Raine, D. J.: 2002, Accretion 
Power in Astrophysics: Third Edition, Cambridge Uni- 
versity Press 

Gaskell, C. M.: 2008, in Revista Mexicana de Astronomia y 
Astrofisica Conference Series, Vol. 32 of Revista Mexicana 
de Astronomia y Astrofisica Conference Series, pp 1-11 

Inada, N., Becker, R. H., Buries, S., Castander, F. J., 
Eisenstein, D., Hall, P. B., Johnston, D. E., Pindor, B., 
Richards, G. T., Schechter, P. L., Sekiguchi, M., White, 
R. L., Brinkmann, J., Frieman, J. A., Kleinman, S. J., 
Krzesinski, J., Long, D. C, Neilsen, Jr., E. H., Newman, 
P. R., Nitta, A., Schneider, D. P., Snedden, S., and York, 

D. G: 2003, AJ 126, 666 

Keeton, C. R.: 2001, Computational Methods for 
Gravitational Lensing, http://xxx.lanl.gov/abs/astro- 
ph/0102340 

Keeton, C. R., Buries, S., Schechter, P. L., and Wambs- 
ganss, J.: 2006, ApJ 639, 1 

Kochanek, C. S.: 2005, in Y. Mellier and G. Meylan (eds.), 
Gravitational Lensing Impact on Cosmology, Vol. 225 of 
lAU Symposium, pp 205-223 

Lawrence, C. R., Elston, R., Januzzi, B. T., and Turner, 

E. L.: 1995, AJ 110, 2570 

Mao, S. and Schneider, P.: 1998, MNRAS 295, 587 
Mathis, J. S.: 1990, ARA&A 28, 37 



6 D. J. E. Floyd, N. F. Bate and R. L. Webster 



He 

A 

m 

D C 


Je 

■ 

■ 

. . .■ ■ 

«, 


Y . 

■ 


z' 


i' 




r' 


g' 

• 




■ ■ ■■ . . !■■ •:■■' ■ 

".'.'■^'■V" ■■-V" ..^ ■ '■'■'•■'.■^^ 

. = 



Figure 1. Magellan IMACS and PANIC imaging of SDSS J0924+0219 taken on 2008 Mareh 21. The anomaly between A and D becomes 
more significant as we observe blucward. The u' image was not used and is included for completeness only. 
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Figure 3. Top: Probability distributions for the radius of the u'- 
band emission region, ctq , in units of the average Einstein Radius 
770, and the power-law index based on our new Magellan data 
only. Contours are shown at 1, 2, and 3a. The dashed lines illus- 
trate four accretion disc models, as described in the main text. 
The small circle indicates the location of the maximum in the 
probability distribution. Bottom: Cumulative probability distri- 
butions for (TO (left) and ( (right). 
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Figure 4. As for Figure [3] but for all available data (this paper, 
Inada et al. 2003, Keeton et al. 2006). The top panel shows the 
probability distributions for u' source radius, ctq and the power- 
law index. The small circle indicates the location of the maximum 
in the probability distribution. The bottom panel shows the cu- 
mulative probability distributions for ag (left) and f (right). The 
alpha-disc prescriptions (f > 4/3) are excluded at the 97% level 
by this combined dataset. 
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Figure 5. Constraints on quasar accretion mechanism in the lit- 
erature (peak probability value ztlcr). Examples shown are: F+09 
(this paper) — SDSS J0924+0219 multi-band single-epoch imag- 
ing; E-l-08 | |Eigenbrod et al.||2008| - QSO 2237-1-0 305 spectro- 
scopic monitoring; E+08ii ( [Eigenbrod et al.||2008[l - a s E-|-08, 
but with no velocity prior; P-l-08 l |Poindexter et al.||2008[l - 
HE 1104-1805 multi-band monitoring; A-l-08 | |Anguita et al.|2008|l 
- QS O 2237-1-0305 multi-band monitoring; B-l-08 | |Bate et al.| 
|2008| l - MG 0414-f 0534 multi-band single-epoch imaging. 



